Search for: All records

Abstract Modifications of the matter power spectrum due to baryonic physics are one of the major theoretical uncertainties in cosmological weak lensing measurements. Developing robust mitigation schemes for this source of systematic uncertainty increases the robustness of cosmological constraints, and may increase their precision if they enable the use of information from smaller scales. Here we explore the performance of two mitigation schemes for baryonic effects in weak lensing cosmic shear: the principal component analysis (PCA) method and the halo-model approach in hmcode. We construct mock tomographic shear power spectra from four hydrodynamical simulations, and run simulated likelihood analyses with cosmolike assuming LSST-like survey statistics. With an angular scale cut of ℓmax < 2000, both methods successfully remove the biases in cosmological parameters due to the various baryonic physics scenarios, with the PCA method causing less degradation in the parameter constraints than hmcode. For a more aggressive ℓmax = 5000, the PCA method performs well for all but one baryonic physics scenario, requiring additional training simulations to account for the extreme baryonic physics scenario of Illustris; hmcode exhibits tensions in the 2D posterior distributions of cosmological parameters due to lack of freedom in describing the power spectrum for $k \gt 10\ h^{-1}\, \mathrm{Mpc}$. We investigate variants of the PCA method and improve the bias mitigation through PCA by accounting for the noise properties in the data via Cholesky decomposition of the covariance matrix. Our improved PCA method allows us to retain more statistical constraining power while effectively mitigating baryonic uncertainties even for a broad range of baryonic physics scenarios. 

ABSTRACT We simulate the scientific performance of the Nancy Grace Roman Space Telescope High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6-yr HLS Reference survey is currently envisioned to image 2000 deg2 in multiple bands to a depth of ∼26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z = 3. The combination of deep, multiband photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g. weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper, we explore multiprobe strategies that can be implemented, given the telescope’s instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the HLS survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the HLS reference survey alone can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets, we assume a flat universe however, and includes realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community-driven effort to simulate and optimize the science return of the Roman Space Telescope. 

ABSTRACT We explore synergies between the Nancy Grace Roman Space Telescope and the Vera Rubin Observatory’s Legacy Survey of Space and Time (LSST). Specifically, we consider scenarios where the currently envisioned survey strategy for the Roman Space Telescope’s High Latitude Survey (HLS reference), i.e. 2000 deg2 in four narrow photometric bands is altered in favour of a strategy of rapid coverage of the LSST area (to full LSST depth) in one band. We find that in only five months, a survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95 per cent of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where the Roman Space Telescope spends 1.5 yr covering the LSST area. For this second scenario, we quantify the constraining power on dark energy equation-of-state parameters from a joint weak lensing and galaxy clustering analysis. Our survey simulations are based on the Roman Space Telescope exposure-time calculator and redshift distributions from the CANDELS catalogue. Our statistical uncertainties account for higher order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modelling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. We find a significant increase in constraining power for the joint LSST + HLS wide survey compared to LSST Y10 (FoMHLSwide = 2.4 FoMLSST) and compared to LSST + HLS (FoMHLSwide = 5.5 FoMHLSref). 

ABSTRACT Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.